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Abstract. We characterized unital quantum channels of single photon polarization 
qubits. The channels are composed of two birefringent crystals and wave-plates, 
where their decoherence properties are controlled. An experimental comparison 
between two different depolarizing configurations was performed using a quantum 
process tomography procedure. The results are with a high fidelity to the theoretical 
predictions. 
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In recent years, within the growing interest in quantum communication, many 
schemes have used the polarization of a single photon to encode a two level quantum 
system - a qubit pQ. Interactions between such a photon and other surrounding systems 
may result in a decoherence process where noise is added to the quantum information 
that is stored in the photon polarization. The possibility to control and characterize the 
properties of this noise enables the study of its effects on various quantum protocols. 
In particular, it will be possible to test the performance of quantum error correction 
schemes and the efficiency of quantum key distribution protocols [21 [3]. The construction 
of a quantum channel that can induce different types of noises is important also for the 
study of a quantum-classical transition that a system may experience (U [5] . 

Implementing a channel with controlled depolarization properties is a challenging 
task, since light interacts very weakly with its environment. Such a channel can be 
realized with the aid of birefringent crystals which entangle the polarization degree of 
freedom with the photon's temporal degrees of freedom [6] or with its spatial degrees 
of freedom [TJ |8]. In this work, we studied a depolarizing channel which couples 
between polarization and time of light of short coherence times [9]. It is composed 
of two birefringent crystals and wave plates. The coupling depends on the relative 
angle between the symmetry axes of the crystals, the wave-plate orientations and the 
initial polarization of the input state. Depolarization occurs because photon detection 
is insensitive to the temporal degrees of freedom. 

The polarization state of a single-photon qubit can be described either by a density 
matrix operator p or, equivalently, by a point in the Poincare sphere. The Cartesian 
coordinates of this point are the Stokes parameters S = {Si, S 2 , S 3 }. Here we use the 
convention that Si represents the linear horizontal and vertical polarizations (\h),\v)), 
5*2 represents the linear polarizations plus and minus 45° (\p) = (\h) + \v))/\/2, \m) = 
(—\h) + \v))/\/2), and S3 represents the circular polarizations (|r) = (\h) -\-i\v))/y2 and 
\l) = (i\h) + \v))/y/2). The length of the Stokes vector D = yj 'Sf + Sf + Sf represents 
the state's degree of polarization. For polarized states, D = 1, while for partially 
polarized states, D < 1. Characterization of the polarization state may be performed by 
several projection measurements with the Quantum State Tomography (QST) procedure 

Consider an arbitrary quantum state p that serves as an input state to a quantum 
channel. The final state p' is derived from the mapping p' = S(p) where S is the 
operation of the channel. The mapping 8 can be uniquely described by the elements of 
the positive and Hermitian matrix \'- 

S(p) = J2XmnE m pEl (1) 

m,n 

where E m are basis elements that span the space of p. The elements of the \ matrix can 
be experimentally determined by a Quantum Process Tomography (QPT) procedure 
[TT] . A QPT procedure for a single-qubit process requires four QST measurements of 
different input states. We note that these four states should not lay on the same plane in 
the Poincare sphere representation. Assuming the process has some special symmetries 
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may reduce the number of QST measurements that are required for the reconstruction 
of the x matrix. A complete positive and trace preserving process is a unital process if 
8(1) — I. A unitality assumption for a single qubit process may reduce the number of 
QST measurements that are needed for the process chracterization to 3. 

The aim of this work is to study theoretically and experimentally the processes 
of two depolarizing channel schemes. The first channel (Scheme I, see figure [1]) is 
composed of two equal birefringent crystals and two half-wave plates. The crystals 
are fixed perpendicularly with respect to each other, and the rotatable wave plates are 
placed before and after the first crystal. The second channel (Scheme II, see figure 
[Q is composed also of two equal perpendicularly fixed birefringent crystals, but here 
a quarter-wave plate is placed in between them. For both depolarizing schemes the 
symmetry axes of the crystals define the linear \h) and \v) polarizations, and the zero 
angle of the wave plates is determined with respect to the horizontal polarization. The 
coupling between polarization and temporal modes of such a channel can be described as 
follows: every crystal induces a temporal delay t between the horizontally and vertically 
components of an input polarized wave packet. This temporal walk-off t depends on 
the crystal length L, its refractive index difference An and the speed of light c such 
that t = L—. In order to achieve a complete separation (i.e. loss of coherence between 
horizontal and vertical wave packets), we require that t > t c , where t c is the coherence 
time of the initial wave packets. We see that if there is no polarization rotation between 
the two perpendicular crystals of our schemes (i.e. a wave plate that is placed in between 
them is set to a zero angle), the second crystal compensates for the time delay that was 
caused by the first one and no depolarization occurs. If there is a polarization rotation 
between the two crystals, two or three different temporal modes will be occupied after 
the second crystal. Tracing out these temporal degrees of freedom while measuring the 
polarization state will result in depolarization. It is easy to calculate the processes that 
these channels induce for any angle of the wave plates because in both schemes the 
polarization is coupled to a maximum of three discrete temporal modes (for details see 
[9]). We emphasize that the wave plates act on any of these modes separately; thus 
they do not couple between the different temporal modes. Theoretical calculations of 
the processes for both schemes have shown that the processes are unital regardless of 
the orientation of any wave plate in the channel configurations. 

Our experimental setup is shown in figure [1] a Ti: Sapphire pulsed laser with 
a 76 MHz repetition rate was frequency doubled in order to generate pulses with a 
wavelength of 390 nm. These pulses were focused into and collinearly down-converted 
in a 1 mm thick type-I BBO crystal. The 780 nm down-converted signal was filtered 
by a dichroic mirror (DM) and collimated with a lens (L2). One photon of the pair 
was probabilistically split by a beam splitter (BS) and was sent to a detector. The 
second photon was directed to the depolarizing channel. For the two depolarizing 
schemes, the birefringent phase of the depolarizing crystals was tuned so that the entire 
depolarizing channel will not affect the state when the wave plate angles are set to zero. 
Characterizing the polarization state of the depolarized photons was performed by wave 
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Figure 1: (Color online) The experimental setup and the depolarizing schemes. Laser pulses of 390 nm 
were focused by a lens (LI) into a BBO crystal and reflected by a dichroic mirror (DM). Down-converted 
signal was collimated by another lens (L2). One photon of the down-converted pair was probabilistically 
split by a beam splitter (BS) and detected by a single-photon detector (DET). The initial polarization 
state of the second photon was determined using half- and qurter-wave plates (A/2, A/4). The signal 
was directed to one of the two depolarizing schemes as shown. QST measurements were carried out 
using wave plates, a polarizer (POL) and another detector. The photons were temporally filtered using 
a 5 nm bandpass interference filters (IF). Spatial filtering was achieved with the coupling of the photons 
to single-mode fibres connected to the detectors. 



plates and a polarizer (POL). Photons were filtered by 5 nm bandpass niters (IF), 
corresponding to a coherence time of t c ~ 180 fs, and then coupled into single-mode 
fibres, leading to single-photon detectors (DET). The detection of the second photon 
was conditioned by the detection of the first one in order to ensure that the depolarized 
signal is truly a quantum state. 

The probabilistic nature of a quantum polarization state, together with systematic 
errors during the projection measurement procedure, may result in an illegal linear 
reconstruction of the measured state or the measured process. A Maximal Likelihood 
(ML) search that restricts the parameters to physically allowed values can be used 
to find the physical representation of p or x that best fits to the tomographic 
data. We reconstructed all of the state density matrices using the ML protocol 
suggested by James et al [12] . As for the % matrices, although ML searches for 
photonic QPT measurements were previously suggested and demonstrated [13j [14], 
we used a simpler ML search that gave sufficient results. According to the Choi- 
Jamiolkowski isomorphism, the linear mapping of a quantum channel represents a 
legal quantum state of a larger dimension [15] . Thus, we can assume that the linearly 
reconstructed process matrix Xunear represents a four-dimensional (4D) quantum state. 
Additionally, we assume that this state was reconstructed from a set of 16 probability 
measurements {pi,P2, ■■iPie} obtained from 16 artificial projection measurements on 
the states |^i6)}- Inverting the linear reconstruction of a density matrix from 

the 16 probability measurements, we calculated these artificial probability values from 
Xunear- Then, we reconstructed a new x matrix using the artificial probability values 
with an ML quantum state search. We emphasis that x might not be the process with 
the best fit to the data, but only a close physical fit to the originally measured data. 

As this ML search requires the translation of the artificial probability measurements 
to a set of artificial counts, there is an extra free multiplying parameter M. For example, 
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Figure 2: (Color online) The first depolarizing scheme: experimentally measured final states and 
QPT in the Poincare sphere representation, (a) Comparison of final state measurements for input 
states of \h) (blue squares), \p) (red triangles), and \r) (green circles) in the range 0° < 9 < 90° and 
the theoretical model. Theoretical curves are presented as solid lines in the range 0° < 9 < 180°. (b-d) 
Mapping of the surface of the Poincare sphere to depolarized wire-mesh ellipsoids that were obtained 
by experimental QPT for the crystal angle values of (b) 9 = 35.3°, (c) 9 = 54.74°, and (d) 9 = 67.5°. 

this parameter can represent the measurement time. We chose M such that the artificial 
set of counts will have the same order of magnitude as the actual qubit measured counts, 
which was about 25,000. Changing M had a minute effect on the x matrix: multiplying 
or dividing M by 10 3 resulted in x matrices whose fidelities with the original x matrix 
(F(XiiX2) = (Tr J 1 ^fx^X2\fX\) 2 ) were higher than 99%. An error estimation of 2% for 
the process matrix fidelities and their eigenvalues resulted mainly from this fact. 

Beginning with the first scheme, we rotated the two half-wave plates in opposite 
directions by an angle of 9/2. This rotation is equivalent to the rotation of the first 
crystal by an angle of 9. We generated initial \h), \p) and |r) polarization states, 
sent them through the depolarizer and performed QST to the depolarized states for 
different 9 angles in the range of < 9 < 90°. The traces of the measured states 
in the Poincare sphere representation are shown in figure [2](a). Solid lines represent 
the theoretical calculations for angles in the range < 9 < 180°. Figures [2]^b)-(d) 
present the measured process of the channel mapping of the Poincare sphere surface 
when 9 = 35.3° = tan -1 (-^), 9 = 54.7° = tan" 1 ^), and 9 = 67.5°, respectively. The 
initial states that were used for the QPT were \h), \p) and |r), where the process of 
9 = 35.3° was reconstruct using the unitality assumption, and the other two processes 
were measured using the additional input state of \v). The first process (figure [2] (b)) 
maps the surface of the Poincare sphere to an ellipsoid with two primary radii of length 
| and the third one has a length of |. The second process (figure [2] (c)) is an isotropic 
depolarizing process that maps the polarized states to a sphere with radius equal to |. 
The third process (figure [2] (d)) corresponds to a mapping of the Poincare sphere surface 
to an ellipsoid with one primary radius with a length of ~ 0.7 and two primary radii 
with a length of ~ 0.15. 

In order to evaluate the channel performance with respect to the theoretical 
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Figure 3: (Color online) The eigenvalues of the \ matrix as a function of the equivalent first-crystal 
rotation angle 9 for the first depolarizing scheme. 



prediction, we compared between the measured % matrices and the theoretically 
calculated Xd matrices. A comparison between the three measured maps that are 
presented in figure [2](b)-(d) and the theoretically calculated maps resulted with fidelities 
F higher than 97 ± 2%. 

A systematic study of the decoherence that the first scheme can induce can be 
performed by exploring the eigenvalues of the 4D \ matrix as a function of the angle 
6. Since they are not affected by rotations and reflections of the input state, these 
eigenvalues may serve as a measure of the pure decoherence types that a quantum state 
may experience while passing through the channel. Using the assumption of unitality, 
we have used the transformed states that are presented in figure |2^a) to reconstruct 
the x matrix for the range of angles < 9 < 90°. The measured eigenvalues of these 
matrices are presented in figure |3] along with their theoretical predictions. We see that 
only three eigenvalues participate in the process, while the fourth eigenvalue remains 
zero for every angle. Two eigenvalues out of these three are always equal. When the 
process is isotropic (8 = 54.7°) the three eigenvalues intersect and are equal to |. 

Moving to the second depolarizing scheme, we characterized the channel process 
for different angles of the quarter- wave plate in the range of < < 90°. The 
QPT was performed by three QST measurements for outputs of the three initial 
mutually unbiased states S a = {^/l/3, 0, —^2/3}, Sb = {-y/l/3, \Jl/2, yl/6}, and 

S c = {Jl/3, — yl/2, Jl/6}. It is worth mentioning that these three states have the 
property that their final output state after passing through the second scheme has 
the same degree of polarization regardless of the angle <fi. Using the assumption of 
unitality, we reconstruct the % matrices of the different processes. The eigenvalues of 
these matrices and their theoretical predictions are presented in figure HI As in the first 
scheme, only three eigenvalues participate in the process, and two out of them are equal 
for any orientation of the quarter-wave plate. Maximal decoherence is obtained when 
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Figure 4: (Color online) The eigenvalues of the x matrix as a function of the quarter-wave plate 
angle <f> for the second depolarizing scheme. 

= 45°, where the eigenvalues are |, |, \ and 0. As can be seen in figure [3j the same 
eigenvalues of the x matrix are obtained also with the first scheme when 9 = 45°. A 
comparison between the measurements of these two processes is shown in figure [5j the 
mapping of the surface of the Poincare sphere according to scheme I process at 9 = 45° 
and the real and the imaginary parts of the x matrix are shown in figure E)^a)-(c), 
respectively. The mapping and the real and imaginary part of the x matrix that are 
shown in figure E)^d)-(f) correspond to the scheme II process at = 45°. Both processes 
map the Poincare sphere into the shape of a disk, but the rotations that accompany the 
processes are different. This is reflected in the difference between the elements of the x 
matrices and in the difference in the orientation of the wire mesh of the mapping. The 
fidelities of these two measured processes to the desired ones are 97 ± 2% and 98 ± 2% 
for the first scheme process at 9 = 45°, and the second scheme process at = 45°, 
respectively. 

In conclusion, we characterized the processes of two depolarizing schemes composed 
of two perpendicularly fixed equal birefringent crystals, and rotatable wave plates. All 
processes were analyzed according to the eigenvalues of their x matrix which serve 
as a measure for the decoherence properties. The schemes realize processes in which 
three eigenvalues of the x matrix differ from zero, and two out of them are equal. A 
comparison of scheme I which contains a half-wave plate in between the crystals and 
scheme II which contains a quarter-wave plate in between them reveals that scheme 

1 offers more possible decoherence types than Scheme II. Both schemes exhibit high 
fidelities to the theoretical predictions. Ongoing research on the decoherence of more 
composite quantum systems is currently in progress using these depolarizing schemes. 

The authors thank Jan Soubusta and Karel Lemr for a helpful discussion, and the 
Israeli Science Foundation for supporting this work under Grant No. 366/06. 
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Figure 5: (Color online) A comparison between the disk mapping processes of the first scheme and 
the second scheme, (a-c) Measured mapping of the surface of the Poincare sphere, real and imaginary 
parts of the x matrix for the first scheme process of 9 = 45°. (d-f) The corresponding results for the 
4> = 45° second scheme process. 
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